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Abstract 

Mathematical arguments are presented that give a unique answer to the question in the title. 
Subsequently, the mathematical analysis is extended using results of detailed model calculations 
that, in addition, throw further light on the consequences of the analysis. Finally, through a 
comparison with various recent studies, many of the latter are given a new interpretation. 
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I. INTRODUCTION 



The response to external electric and magnetic fields provides a fundamental tool for 
studying and altering the properties of materials with numerous attendant applications. In 
particular higher-order responses allow for 'manipulating light with light.' Thus, there is 
considerable interest in identifying molecular systems with large non-linear responses. 

One approach in this direction is based on push-pull systems, i.e., chain-like molecules 
with an electron-donor group at one end and an electron-acceptor group at the other (see 
Fig. [1]). When the backbone is a ^-conjugated oligomer the 7r electrons of the backbone may 
respond easily to perturbations like those of the substituents and/or external fields. Due 
to the donor and acceptor groups a large electron transfer, and, accordingly, a large dipole 
moment can occur and one may hope for large responses of the dipole moment to external 
fields. For these 7r-conjugated systems, each circle in Fig. [T]could be, for example, a vinylene 
group, a phenylene group, a methinimine group, or combinations of those. 

If the push-pull system is sufficiently large, we may split it into three parts, i.e., a left 
(L), a central (C), and a right (R) part as shown in Fig. [TJ Electrons of the central part 
are assumed to be so far from the terminations that they do not feel the latter (or, more 
precisely, the effects of the terminations are exponentially decaying in the central part). 

The dipole moment, jl, is useful in quantifying the response of the system to an external 
electric field, 



Here, E m (u s ) is the mth component (i.e., x, y, or z) of the external field with the frequency 
uj s and u is the frequency of the response of the molecule to the field. The u n summations 
go over all the frequencies of the applied field. /j,\°\u) is the dipole moment in the absence 
of the field which vanishes for u 7^ 0. Moreover, ctij(u)]±wx) is the linear polarizability, 
and Pijk(w, ±wi, ±^2), lijki(w, i^i, i^2j ±^3)) ••• are the first, second, ... hyperpolariz- 
ability. Sum rules require that these quantities can be non-zero only if the frequency of 
the response, u, equals the sum of the frequencies (eventually multiplied by —1), i.e., for 




+ eH H lijki(u; ±uji, ±u 2 , ±uj 3 ) ■ E j (Lj 1 )E k (u 2 )Ei(u; 3 ) H . 



(1) 
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Jijkliu; ±^i, ±^2, ±^3) we require uj = ±o?i ± co>2 ± cj 3 . 

In the present paper we focus on static external fields, in which case uJi = 0. Furthermore, 
we shall study a neutral system, although our arguments also are valid for charged systems as 
long as the extra charge is localized to the terminations. We let p(r) be the (field-dependent) 
total charge density (i.e., the sum of the nuclear and electronic charge densities), and choose 
the long axis to be z. Then the component of the total dipole moment that is of interest 
here, namely z, is given by (omitting its argument, uj) 

p z = J p(r)zdf= J p(r)zdf + J p(r)zdf+ J p(r)zdr, (2) 

where we have split the integral into contributions from the left, central, and right regions 
of the chain. The central region consists of identical neutral units. We can, therefore, write 



p(r)zdf=K c pc, (3) 

c 

where Kc is the number of units in C and pc is the z component of the dipole moment 
of one of these units. In order to evaluate the other two contributions to the total dipole 
moment in Eq. (|2J) we define a 'typical' center for each term, i.e., Rr and Rl (these could, 
e.g., be the center of mass of the right and left parts, respectively), and let Zr and Zl be 
the z components of these vectors. Since the chain is neutral we, then, obtain 



p(r)zdr+ I p(r)zdr = (Zr — Zl) / p(r)dr+ / p(f)(z — Z£)dr-\- / p(r)(z — ZjAdr. (4) 

L JR JR JL JR 

The first term on the right hand side describes the contribution to the dipole moment 
associated with electron transfer from one end to the other. This term grows linearly with 
chain length (due to Zr — Zl) as does the term in Eq. On the other hand, the last two 
terms in Eq. (jlj) describe local dipole moments that arise from the electron distributions 
within the two terminal regions and they are independent of the chain length. 

This discussion suggests that donor /acceptor (=D/A) substitution at the ends of long 
chains may change the charge distribution in R and L so as to strongly enhance the dipole 
moment and, consequently, produce a particularly large change in the dipole moment when 
the system is exposed to an external electric field. Therefore, very many stjidjes i Jmw i j3een 
devoted to push-pul . 



systems function of increasing length (see, e.g., 

QQQQIii 15 3 3H ) . 

Not only the electrons but also the structure (phonons) will respond to a static electric 
field. We will demonstrate that, for sufficiently long chains, the electronic response per unit 



of a push-pull system (with structural relaxation taken into account) becomes independent 
of the donor and acceptor groups, implying that the materials properties cannot be improved 
upon substitution. Our mathematical arguments for this finding are presented in the next 
section, and in Sec. IHII we illustrate and analyse the results through calculations on a model 
system. The particular case of inversion symmetry is discussed in Sec. IIVI where we also 
make a comparison with previous results. Finally, a summary is provided in Sec. |Vj 

The arguments we present are related to those originally given by Vanderbilt and King- 
Smith for an extended system in the absence of an external field. They argued that the 
permanent polarization (i.e. dipole moment per unit length) is a bulk property. 1 ^ Very re- 
cently, Kudin et alM proved that the permanent polarization is quantized for D/A substi- 
tuted systems. Neither of these works considered the induced polarization or the structural 
relaxation due to an external field. Finally, in a recent paper we presented some of the 
arguments behind the present work but did not analyze the predictions as we do here using 
a model system.— 

II. CHANGES IN THE CHARGE DISTRIBUTION UPON SUBSTITUTION 

By replacing some (groups of) atoms with others at the chain ends, the electronic orbitals 
with components near the ends will change. Since the set of electronic orbitals is orthonor- 
mal, all other orbitals will change as well. Accordingly, the charge distribution may change 
everywhere due to the substitutions. 

When an electrostatic field is applied as well, each orbital will respond to the field. Since 
the orbitals will have changed due to the substitution, so will their responses to the field. 
Furthermore, the structural responses due to the field will also depend on the substitution 
at the ends. Therefore, the dipole moment can depend upon both the substitution and the 
field. From these arguments there is no reason to believe that /i (0) /A, a/N, (3/N, j/N, 
. . . (with N being the number of repeated units) will be independent of the substitution. 
However, we shall argue here that the charge 



in Eq. (j3J) can change, at most, by an integral number of elementary units for different D/A 
substitutions at fixed external static field. Our proof is a generalization of arguments due 




(5) 
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to Vanderbilt and King-Smith 19 (see also 20]), and was previously proposed by the present 
authors.— It will be verified here by calculations on a model system and given a thorough 
analysis on that basis. 

For a given system (with specified geometry), and value of the external field, we transform 
the set of occupied orbitals into a set of orthonormal, localized functions. Those functions 
ascribed to C will be similar to the Wannier functions of the infinite periodic system. The 
localized orbitals will be centered in one region, but may have tails sticking into another 
region. We assume that the terminal regions are large enough so that any functions centered 
therein, which differ from those of C, are exponentially vanishing in C. On the other hand, 
those functions ascribed to C, but centered on units next to L or R, will likely have tails 
extending into those regions. 

The density matrix can then be written in block-diagonal form with three blocks, one for 
each of the three regions. Since the density matrix is idempotent, each block will be so, too, 
and there will be an integral number of electrons associated with each of the three sets of 
functions. That is to say, the number of electrons associated with the functions centered in 
the two end regions is integral. Accordingly, any non-integral part of q is associated with 
the tails of the functions in C that extend into R, which, per construction, is independent 
of the terminations, i.e., also of D/A substitution. 

We conclude that, for different terminations, q can change only by an integer. This is 
valid for long chains and all fields. Therefore, the electronic response per unit of the chains 
to the field, with or without nuclear response, is independent of termination. The only 
possible change for different terminations is that q may jump by an integer for different field 
strengths. In fact, our numerical studies on a Hiickel-type model will confirm this prediction. 
Of course, in ab initio calculations, there may also be a jump due to changing the basis set 
or the method (e.g. Hartree-Fock vs. Kohn-Sham DFT). 



III. ILLUSTRATING AND ANALYZING THE RESULT 

In order to explore in detail the predictions from above we studied a Hiickel like model 
for long, finite (AB)2K+i chains. In our model, we use a basis set of orthonormal atomic 
orbitals (AOs) with one AO per atom. The system has one electron per atom, and the nuclei 
are given charges of +1 whereas the electronic charge is set equal to —1. (All quantities 
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are expressed in atomic units in this paper.) Given that Xn is the AO of the nth atom 
(n = 1,2, ...,AK + 2) and h is the Kohn-Sham or Fock single-electron hamiltonian we 
assume that only (Xj\h\Xj)i (Xj\h\Xj±i) i an d (Xj\h\Xj±2) are non-vanishing with values 

(X2 P +i|%2 P +i> = e 

(X2p\h\X2 P ) = -e 

(Xj\h\Xj+i) = ~[ti ~ ai(z j+ i - zj)] 

(Xj\h\Xj+2) = -[t 2 - a 2 {z j+2 - Zj )}. (6) 

Here Zj is the position of the jth atom. Different donor and acceptor groups are modeled 
by modifying the on-site energies of the terminating atoms and/or the terminating hopping 
integrals, 

(Xi\h\Xi) = eo + £L 
(X4K+2\h\X4K+2) = -e + e R 

(Xl|^IX2> = - (*1(Z2 ~ Z X )} +t L 

(X4K+l\h\X4K+2) = ~[tl- «l(^4X+2 ~ Z 4K+1 )] + t R . (7) 

Finally, we assume that 

(Xj\z\Xk) = 6j,kZj- ( 8 ) 
In order to analyse the results we, first, define a reference structure for which the position 
of the nth atom is 

4 0) = \ (n - 2K - - (-1)^ . (9) 

Here a is the length of the unit cell for an infinite, periodic system with the same electronic 
interactions and no external field. Subsequently, we define for each atom 

u n = z n - | (n - 2K - 

Az n = Zn-Z®. (10) 

The total energy is written as the sum over occupied orbital energies (multiplied by 2 due 
to spin degeneracy) augmented by a harmonic term in the nearest- and the next-nearest- 
neighbour bond lengths, 

occ ^ 4K+1 AK 4K+2 

E tot = €i + y ( Z P+1 - Z vf + y H(^P+2 - Z p f - E DC Z p . (11) 

i=l p=l p=l p=l 



i?DC is the strength of the electrostatic field. For the infinite, periodic chain without an 
external field, the lowest total energy corresponds to a certain lattice constant a and 

u n = (-l) n+ V (12) 

The force constants k\ and k<i are determined so that a and uo take certain chosen values. 
With 

4K+2 

A = Yl C niXn (13) 
n=l 

being the ith orbital (ordered according to increasing orbital energy) we calculate the Mul- 
liken charge on the nth atom for field E DC as 

2K+1 

q n (E DC ) = 1-2^ \C m \ 2 (14) 



i=l 



which leads to the dipole moment 



4K+2 

Hz= z n q n (E DC )- (15) 



n=l 

The charge transfer is given through 



4K+2 

q = Yl Qn(E DC )- (16) 

n=2K+2 

We also define 

Aig n (£ D c) = q n (E DC ) - q n (0) 

A 2 q n (E BC ) = q n (E DC )-qn(0). (17) 

where q n (0) is the charge for the infinite, periodic chain in the absence of the field. A 2 g„(i?Dc) 
quantifies the effects on the charge distribution of the push-pull chain due to including the 
field, whereas Axg n (_E DC ) includes effects both from the field and from the terminations. 
Note that Aiq n (Ei)c) — A 2 g n (-EDc) gives the field-independent effect of the terminations. 
Finally, it turns out to be useful to define the center and width of the ith orbital according 
to 

4K+2 



-nil 



AC, 



n=l 

4K+2 "I 1 / 2 

y ] ( z n Ci) |C m | 
n=l 



|2 



(18) 
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which is consistent with Eq. (jSJ). 

We performed calculations for six different terminations specified by (e^, €r, tL, £r). The 
results are summarized in Figs. [2J El HI El and O Since our model is that of a finite chain 
with two different types of atoms, A and B, the Mulliken charges in the central region take 
two values. This is clearly recognized in the presentation of q n in Fig. [2] for E^c — —0.015. 
In Fig. [2] it is also seen that near the ends, the Mulliken charges differ from the values of the 
inner part and, moreover, these charges depend sensitively on the terminations. For the field 
strength E-^q = —0.015 these findings are only marginally modified compared to those of a 
vanishing field (not shown). From Aiq n we see that the combination of electrostatic field 
and termination leads to an internal polarization of each unit in C. Actually, Aiq n shows a 
reduced internal polarization compared to A 2 q n - Thus, terminating the chain reduces the 
effect of the field in that regard. Whereas A2q n contains information about the field-induced 
charge redistributions, Aiq n contains additional information about the (field-dependent) 
effects of the terminations. For E^c = —0.015 the field-induced charge redistributions are 
smaller near the terminations than in the central parts. 

For the larger field, Edq = —0.03, in Fig.[3]the identification of the central region becomes 
much more difficult and, as we shall see below, electrons are transferred from one end to the 
other. Moreover, in this case the field perturbs the system so strongly that the effects of the 
field are stronger than those of the terminations. This can be seen from the fact that Aiq n 
and A2q n are very similar. 

The structure also depends upon the termination. For the intermediate field of Fig. [2] 
(and for zero field as well) the atomic coordinate u n is nearly constant in C but varies 
considerably near the ends where its value depends on the termination, as was the case for 
the atomic charges. For the higher field in Fig. [3] it appears as if no central region can 
be identified from this parameter. However, the fact that Az n is essentially linear for the 
innermost atoms implies that there is a well-defined, repeating structure in C with a lattice 
constant differing from that of the field-free case. 

Fig. H] shows that the charge transfer, q, is independent of termination (though not inde- 
pendent of the field), with the exception of jumps by (even) integers. (The integers are even 
because we have not allowed for spin polarization.) However, the charge distribution inside 
R or L does depend on the terminations and, as a consequence, the dipole moment does as 
well. On the other hand, the variation of fi z /N as a function of E^c for different termina- 
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tions follows parallel curves, implying that the (hyper) polarizabilities are independent of the 
terminations. In fact, a least squares fit yields the values (including maximum deviations): 
Ho/N = 0.3245 ± 0.0023, a/N = 1.677 ± 0.013, (3/N = 18.17 ± 0.20, and 7/iV = 606.9 ± 9.1 
for all six terminations. 

As a function of field fi z is discontinuous and the power series expansion is valid only 
up to the field where the discontinuity occurs. Once such a discontinuity has been passed, 
the dipole moment depends more strongly on the field. This means that the only way of 
increasing the responses of long push-pull systems to DC fields is to design chains for which 
the integral electron transfers occur at low fields. 

At a given field the size of the chain for which jumps in the charge q (i.e. Zener tunneling) 
take place depends on the terminations (cf. Fig. In the shortest chains, for which Zener 
tunneling does not occur, \x z follows parallel curves as a function of chain length, N = 2K+1, 
for different terminations. This means that the dipole moment and (hyper)polarizabilities 
per unit become independent of termination. However, as seen in Fig. [5J the slope of these 
curves increase after Zener tunneling has taken place, implying that the dipole moment 
increases. Assuming that the field- dependence of the dipole moment likewise increases, this 
suggests that the polarizability and/or hyperpolarizabilities per unit may increase for D/A 
substituted systems after an integral number of electrons has been transferred from one end 
to the other. 

In Fig. [6] we show an example of what happens to the molecular orbitals when the jumps 
take place. Calculations were performed for field strengths between —0.0340 and —0.0485 in 
steps of —0.0005, but in the figure we only show the results for fields where Zener tunneling 
occurs. In all cases, the curves vary smoothly as a function of field strength. At the lowest 
two fields, the occupied orbitals closest to the Fermi level have a center in the left part 
(Ci < 0), whereas the unoccupied orbitals closest to the Fermi level are centered in the right 
part. At the field E-qq ~ —0.0375, two electrons (one per spin direction) are transferred 
from one side to the other, which again happens at a larger field (-Edc — —0.0475). In the 
first case, we observe the occurrence of two new, very localized, orbitals close to (but not at) 
the Fermi level. The energetically lower (i.e occupied) one is localized towards to the chain 
end on the right side while the other (unoccupied) is localized towards the chain end on the 
left side. Accompanying this interchange is a similar interchange of two rather delocalized 
orbitals, both of which are further away from the Fermi level and centered closer to the 



9 



middle of the chain. Again, at the second electron transfer a pair of new, rather localized, 
orbitals near (even closer to) the Fermi level show up towards the chain ends, and also this 
transfer is accompanied by some reorganization of the other orbitals. Finally, Fig. O also 
shows an example of a reorganization of the orbitals, i.e., for a field around Edq = —0.0430. 
Here, one localized, occupied orbital interchanges order with an adjacent (in energy) more 
delocalized orbital, but otherwise no further significant changes are observed. 

IV. INVERSION SYMMETRY AND COMPARISON WITH PREVIOUS RE- 
SULTS 

Before proceeding to compare with previous results we develop an interesting consequence 
of our findings with regard to inversion symmetry. The same arguments can be applied for 
a system containing a mirror plane perpendicular to the chain axis, but here we shall for 
the sake of simplicity restrict ourselves to the case of inversion symmetry. Suppose the long 
oligomer of interest contains a central region made up of units with inversion symmetry. 
Even if the central part does not have inversion symmetry, it may be possible to create 
such with the addition of appropriate terminating groups. This is, for example, the case for 
oligomers of thienyleneethynylenes and thienylenevinylenes that were studied by Geisler et 
al— Many of the systems of interest fall into one of these two categories. Since, according to 
our findings, D/A substitution cannot affect the (hyper)polarizabilities per unit, the latter 
must vanish even if the symmetry is not preserved. For instance, modifying the terminations 
of the systems of Geisler et al. so that inversion symmetry no longer exists cannot result in 
a non- vanishing fl/N if the chains are sufficiently long. 

A large fraction of previous observations are for systems of the type described in the 
preceding paragraph. Some of these cases are discussed below along with others pertinent 
to our findings herein. We now briefly consider, in particular, the works mentioned in the 
Introduction. 

In their combined experimental and theoretical study on some push-pull oligoenes, Mey- 
ers et air- observed a 'negligible charge transfer all the way from the donor to the acceptor', 
which implies that q is independent of the termination. On the other hand, in their theo- 
retical study Tsunekawa and Yamaguchi 2 examined shorter, nitrogen-containing push-pull 
oligomers. They noted that these systems are interesting from the perspective of maximizing 
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(3, but our results establish that, for such to be true, the systems must be short enough so 
that our approach is inapplicable. This serves to highlight the point that apparent, but not 
real, discrepancies can occur due to shortness of the chain length. 

Marder et al.— presented an approach for unifying the description of linear and nonlinear 
polarization in organic polymethine dyes. It has since been shown that their analysis is 
invalid if phonons are taken into account.— Here, however, we emphasize that the conclusions 
they draw regarding (3 can, again, hold only for systems that are too short for our treatment 
to apply. 

Clearly, the chain length required for validity of the treatment given here is an important 
issue. In Fig. the dipole moment is converged for chains with some 20 units. However, 
this may be an artifact of our simple Hiickel model. In an experimental study 4 and in 
several computational studies,- 1 *^ 1 ^ 4 - 1 ^ the second hyperpolarizability per unit was found to 
converge considerably slower which, in fact, agrees with our own earlier findings.— Thus, 
when focusing on higher-order non-linear responses quite large chains may be required for 
the results of the present work to be relevant. In shorter push-pull systems (for instance 
those considered by Geisler et a/.- 1 ^ or by Morley et al.Z&) D/A substitution can have an 
influence on the response. 

As shown numerically by Champagne et al.,— (3/N also converges relatively slowly as 
a function of size. They considered D/A substituted oligomers of polymetheimine [also 
called polycarbonitrile, (CHN)J. This system has a zigzag backbone of alternating C and N 
atoms with alternating bond lengths. Without the bond length alternation it would, at least 
hypothetically, be possible to choose donor and acceptor groups so that the overall system 
is centrosymmetric. Even if chemical arguments imply that this structure is unrealistic, a 
non-zero value of f3/N for long chains should be ascribed, strictly speaking, to the bond 
length alternation. 

Polyphenylenes and polypyridines have been studied by Zhang and Lu.— They focused 
on a and 7 as a function of the length of a closed ring for each system and applied a finite- 
field approach in their calculations. Unfortunately, as we have shown earlier (see, e.g., j^H), 
this approach will never converge to the results for the infinite, periodic chain. Nevertheless, 
although (3/N will vanish for the polyphenylenes, we predict that a non-zero value will occur 
for both short and long oligomers of the polypyridines. 

For the D/A substituted polyenes studied by Champagne et al.— our analysis confirms 
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their findings, i.e., that j3/N will vanish for sufficiently large chains. Their numerical results 
indicate that f3/N goes through a maximum and that convergence to the infinite chain result 
for larger N is slow. 

Even the polarizability, a/N, and the permanent dipole moment, fx^/N, may converge 
more slowly, as a function of chain length, than predicted by our simple model. This is, for 
example, the case for the systems investigated by Smith et alM and by Kudin et alM 

In a recent study, Botek et al— compared finite oligomers of [iVjhelicenes and 
[iV]phenylenes that possess a helical structure for N larger than roughly 6. By making ex- 
plicit use of the helical symmetry of the central region we predict that, when those systems 
are sufficiently long, D/A substitution will not be able to modify the electronic responses 
to static fields. The fact that Botek et al. find changes upon D/A substitution implies that 
the chains of their study are not converged to the long chain limit. 

V. SUMMARY 

As long as the applied field is not so strong that an integral number of electrons is 
transferred from one end to the other, the answer to the question of the title is clearly: there 
can be no change. This comes from our mathematical analysis in Sec. [TTJ which generalizes 
treatments presented previously by Vanderbilt and King-Smith^ and by Kudin et al.,— who 
considered only electronic polarization in the absence of an external electrostatic field. It is 
also in agreement with our own earlier prediction. 21 

Calculations on a model system confirm the basic result and shed light on the nature of 
the end-to-end charge transfer. Although the end charges, permanent dipole moment, and 
structure depend sensitively on the terminations neither the amount of charge transferred 
nor the (hyper)polarizabilities per unit do so. The field and/or chain length at which the 
charge jumps take place also depend on the terminations. Each jump is associated with 
an interchange of occupied and unoccupied molecular orbitals that are well-localized in the 
chain end region. These orbitals are close to but not at the Fermi level. There is also an 
accompanying orbital reorganization. 

One consequence of our finding is that long unsubstituted chains which have inversion 
or mirror symmetry, or can be made symmetric by substitution, must have a vanishing first 
hyperpolarizability per unit. Experimental and theoretical determinations are consistent 
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with this fact, although apparent contradictions can occur for short chains. 
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FIG. 1: Schematic representation of a long finite chain separated into a central region and two 
terminal regions. The lower chain shows what happens when a donor and an acceptor group have 
been added at the terminations. Each dot represents a group of atoms. 
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FIG. 2: Results from the model calculations for chains with K = 40, i.e., 162 atoms and an external 
field of -0.015 a.u. The first three rows show quantities related to the charge distribution, i.e., those 
of Eqs. (fT3|) and (fTT|) . whereas the two lowest rows show quantities related to the structure, i.e., 
those of Eq. (fT0|) . In all cases, eo = 0.5, t\ = 2.5, a.\ = 4.5, ti = 0.5, and = 0.2. Moreover, k\ 
and &2 were chosen so that the optimized structure for vanishing field and an infinite periodic chain 
has a = 2.1 and uq = 0.05. Open circles, closed circles, crosses, open squares, closed squares, and 
open triangles show results for which (e^, e#, ix,, tn) has been set equal to (0,0,0,0), (0,0.5,0,0), 
(0.5,0,0,0), (0,0,0,0.8), (0,0,0.8,0), and (0.5,0.5,0.8,0.8), respectively. 
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FIG. 3: Results from the model calculations for chains with K = 40, i.e., 162 atoms and an external 
field of -0.03 a.u. The left and middle columns show results for the two ends whereas the right 
column shows the results for the whole chain. The presentation is as in Fig. [2j 
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FIG. 4: The charge transfer, q, and the dipole moment per unit, fi/N with N = 2K + 1, for the 
same systems as in Fig. [2j The right panels show a magnification of the left ones. 
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FIG. 5: The charge transfer, q, and the dipole moment \x with N = 2K + 1, for the same systems 
as in Fig. [2] as functions of N. The two lower rows show a magnification of the two upper ones. 
The field strength has been set equal to 0.03, 0, and —0.03 for the left, middle, and right columns, 
respectively. 
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FIG. 6: Q an d A£j as a function of orbital index for a given push-pull chain with different fields. 
The vertical dashed lines separate occupied and unoccupied orbitals; K = 20; and the parameter 
values are the same as in Fig. [2] except that only the first termination is considered. The field 
strength is given in each panel. Since A£j > the two curves are readily distinguished. Moreover, 
the panels for -Edc = —0.0340 and —0.0485 have the same orbital index scales, as is the case for 
the panels for £ D c = -0.0370, -0.0375, -0.0425, -0.0430, and -0.0435, and for the panels for 
E DC = -0.0470 and -0.0475. 
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